Two-dimensional small-world networks: navigation with local information 
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Navigation process is studied on a variant of the Watts-Strogatz small world network model 
embedded on a square lattice. With probability p, each vertex sends out a long range link, and the 
probability of the other end of this link falling on a vertex at lattice distance r away decays as r~ a . 
fvj ■ Vertices on the network have knowledge of only their nearest neighbors. In a navigation process, 

messages are forwarded to a designated target. For a < 3 and a/2,a scaling relation is found 
between the average actual path length and pL, where L is the average length of the additional long 
range links. Given pL > 1, dynamic small world effect is observed, and the behavior of the scaling 
function at large enough pL is obtained. At a = 2 and 3, this kind of scaling breaks down, and 
different functions of the average actual path length are obtained. For a > 3, the average actual 
path length is nearly linear with network size. 
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PACS numbers: 89.75.Hc, 84.35. +i, 87.23. Ge, 89.20.Hh 



I. INTRODUCTION 



The famous experiment of messages being forwarded to a target among a group of people, carried out by Milgram 
[1] in the 1960s, and also by Dodds et al. in 2003 [2] on a larger scale, reveals the existence of short paths between pairs 
of distant vertices in networks that appear to be regular (i.e. the small- world effect [1,3-6]). One of the important 
quantities that characterize this small-world effect is the average shortest path length (d) between two vertices. On 
small world networks, this value grows very slowly (relative to the case of a fully regular network) with the network 

O , size N. Recent empirical research has shown that a great variety of natural and artificial networks with their structure 
dominated by regularity are actually small worlds, and their average path lengths grow as In TV [7], or more slowly. 
T — l | (See Refs. [8] for more reviews.) 

^ ■ An alternative issue revealed by experiments, but less obvious, is about the realistic process of passing information 

on small world networks. This kind of information navigation has been studied by Kleinberg [9]. This process goes 
dynamically: When a message is to be sent to a designated target, each individual forwards the message to one of 
its nearest neighbors (connected either by a regular link or a shortcut) based on its limited information. Without 
information of the whole network structure, this actual path is usually longer than the shortest one given by the 
topological structure. While (d) is the average shortest path length, the average actual path length (I) is the average 
number of steps required to pass messages between randomly chosen vertex pairs. As (d) is usually referred to as the 
diameter of the system, in the rest of this paper (I) shall be taken as the effective diameter, and (d) < (I). It has been 
noted that the topology of the network may significantly affect the behavior of (I) [9,4,10-12]. In other words, it may 
determine the efficiency of passing information. 

Based on Milgram's experiment, Kleinberg studied navigation process on a variant of the Watts-Strogatz (W-S) 
small- world model [3] on&n NxN open regular square lattice. Each vertex sends out a long range link with probability 
p, and the probability of the other end falling on a vertex at Euclidean distance r e away decays as r~ a . Kleinberg 

$-£ ' studied (I) when each vertex sends out one long range link, and proved a lower bound of (l P =i) = cN@( a \ When 
a = 0, the long range links are uniform, and (Z Q =o,p=i) oc N 2 / 3 was obtained, de Moura et al. [11] studied (l) on the 
D-dimensional W-S model with a = and varying p, and obtained (Z a= o,p=i) oc N 1 ^ D+1 \ and thus (Z Q =o,p=i) oc N 1 / 3 
in the two-dimensional case. In the more recent work of Zhu et al. [10] on the one-dimensional case, the variance 



of (l a ,p) with both a and p was studied, and scaling relations were shown to exist. For the studies of the searching 
processes on other different networks, see Refs. [4,13]. 

In this paper, we systematically investigate the navigation process on a two-dimensional variant of the W-S network 
model [3,14]. We study the behavior of (I) by first working out the scaling relations in the two-dimensional case. Our 
result also provides new understanding of the scaling analysis in Ref. [10]. In Sec. II, the model used here is 
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constructed and the navigation process is described, and then the average actual path length (I) is obtained with 
some approximation based on a rigorous treatment. Following that, in Sec. Ill, the dependence of (l) on N,a, and p 
are presented based on scaling relations. Special attention is paid to the cases studied in the works of Klcinbcrg [9] 
and de Moura [11]. Summary and discussions can be found in Sec. IV. 

II. THE CONSTRUCTION OF THE NAVIGATION MODEL 

Our model starts from a N x N two-dimensional square lattice. With periodic boundary condition, the lattice 
distance between two vertices (x, y) and (x 1 , y') can be written in a two-dimensional fashion as 

*•(*,„),(*',„') = Az + Ay, Ax = N/2-\\x-x'\-N/2\, Ay = N/2 - \\y - y/\ - N/2\ . (1) 

This value is actually the length of the shortest path connecting these two vertices through only regular links. To 
generate a small world, with probability p (0 < p < 1) each vertex sends out an additional link to another vertex 
(excluding its original nearest neighbors) 1 . If this other vertex is selected at random, then we are creating a small- 
world network with random shortcuts. Based on realistic considerations (for example, people tend to be brought 
together by similar interest, occupation, etc.), we shall add the shortcuts in a biased manner [9,14]: If the shortcut 
starts from vertex i(Xi,yt), the probability that vertex j{xj,yj) is selected as the end depends on the lattice distance 
between them in the following way, 



" \ r ( x ,y),(. x ',y')) ^ r (x,y),(x',y') J 



where a is a positive exponent and 



(x" ,y")=£(x,y),{x±l,y),(x,y±l) 



(x,y),(x",y") 



is the normalization factor. 

In the model described above, the navigation process can be simulated with the so-called "greedy" algorithm [9]: 
Without loss of generality, suppose the target is vertex (0,0). At each step, the current message holder, vertex 
n(x n ,y n ), passes the message through one of its regular or long-range links. Based on its limited local information, 
this link is believed to bring the message the closest to the target. Based on this algorithm, the actual path length 
(l(x n ,y n )) can be obtained after taking an ensemble average over all possible realizations of the network (with a set 
of fixed parameters, p, N, a, etc.). 

In the simplest case, we suppose that each vertex has information of only the vertices that can be reached within 
one step, and do the calculation as the following: (1) If the current message holder is vertex (0, 1), we simply have 

<Z(0,1)} = 1. (2) 

It is the same for the other nearest neighbors of the target (0,-1), (1,0), and (—1,0). (2) There are 8 nodes with 
lattice distance 2 from the target: (0, ±2), (±2,0), (±1,±1), and (=pl,±l). If the current message holder is, for 
example, (0,2), then with probability 

W(o,2)->(o,o) = 1 - f 1 -"P— 

it is directly linked to the target via one shortcut, which means the message is sent directly to the target with this 
probability. On the other hand, the probability that the message is forwarded along a regular bond is 



W reg — 1 - ^(0,2)^(0,0)- 



Thus 



*It should be noted, though, that we do not consider direction of links in the following discussions. 



(I (0, 2)) = 1 x W (0 ,2)^(o,o) + 2 x W reg . (3) 

The calculation is the same for the other 7 nodes mentioned above. (3) In a general case, the message is held by 
vertex (x, y). W/ xy \^/ x i y >\ denotes the probability that the message is forwarded in the next step to a vertex (x 1 , y'), 
which must be nearer to the target than (x, y) by at least lattice distance 2. If the message holder is not able to find 
a shortcut, the message will be forwarded along a regular link with probability 

W reg = l- J2 W (x,vWx',v')- ( 4 ) 

(x',y') 

For example, if the vertex (3,2) passes the message through a regular link, it will randomly choose (2,2) or (3,1), 
which in the following will be denoted by (x reg ,y reg ). 
Now, with this set of probabilities W, we obtain, 

(!(a;,|/))=W r ( Ii!/ ^(o,o)+ J^ W {x ^ y) ^ (x ^ y/) [1 + (l(x',y'))} +W reg [l + (I (x reg ,y reg )}} . (5) 

(x',y') 

Considering that p/A is a relatively small quantity, W^a:,^)— >(x',3/') can be expressed as [10] 

r -a 

W — o- (x,y),{x' ,y') , . 

W( x ,y)^{x',y') ~ IP J , (D) 

where we have used the fact that r(x,y),(x',y') = r (x'.y').(x,y)- Then it is an easy task to obtain W reg from Eq. (4). 

Recall that we have the definition of the average shortest path length (d) = n 2 (n 2 -d X^i dij, where d^ is the 
length of the shortest path between vertices i and j. By contrast, the average actual path length can be defined as 

Further, with vertex (0,0) being the target, we group the other nodes according to their lattice distance from the 
target, and in the following we shall also discuss the function (I (n)), which for each value of n is obtained by averaging 
all nodes (x,y) with r/j^mm = n. 

III. FEATURES OF THE NAVIGATION PROCESS 

The average actual path length depends on multiple parameters. Here we take into consideration varying N, a, 
and p, but keep the range of view of each vertex limited to its nearest neighbors. Our discussion of the navigation 
process starts from looking for the basic scaling relations. 

Scaling relation is not new in the theories of small world effect. Actually, it plays a central role in the current 
theoretical framework. In 1999, Newman [15] showed that in the W-S model with uniform shortcuts the average 
shortest path length is a function oipN, and it sharply decreases when pN becomes larger than 1. Newman noticed 
that pN is simply the expected number M of long-range links. The threshold of small world behavior is M = pN > 1, 
that means the network becomes a small world when there is more than one long-range link. When the model network 
is generalized, this interpretation shall be generalized as well. For example, in a discussion of the scaling relations 
in the problem of dynamic navigation, Zhu et al. [10] considered inhomogeneous long range links (the probability of 
linking two nodes falls when their lattice distance increases). In their study, the dynamic small world behavior is 
switched on when ML' > 1, where M is the number of long-range links and L' is the average reduced link length (the 
average length of long range links L divided by the system size) . Although they focused on different aspects (static 
and dynamic) of the small world effect, we can still compare these two versions of scaling relations. Because in the 
model studied by Newman, L' ~ 1/4, it is consistent with the interpretation of Zhu et al. Actually, as we shall see 
below, the interpretation of Zhu et al. can be developed as well, when a more general model is considered. 

In the Introduction we have defined (l) as the effective diameter, in the following we will use (I') = (l)/N as the 
reduced effective diameter. If the network is regular, (/') will appear as a constant. Bearing this in mind, we first look 
at the results shown in Fig.l. For each value of a (with exception at a — 2, as will be shown below), (/') appears as a 
function of pL, where L is the average length of long-range links. Thus our study clearly supports an interpretation 
different from that in Ref. [10]: Instead of ML', the parameter should be pL. In the one-dimensional case, this equals 



ML', and thus consistent with Rcf. [10]. When pL <C 1, (l 1 ) — > 0.5, and (Z) a TV, indicating that the network is 
virtually regular, and when pL increases beyond 1, the system begins to show a dynamic small world behavior. 
However, we find interesting exceptions at a = 2 and a = 3. As shown in Fig. 2, at a = 2, exp ( (I') x pL J is a linear 

function of pL/\nN for pi significantly larger than 1. Due to this extra factor of 1/lnJV, there is no way that (/') 
can be written as a function of pL for a = 2. Because at a = 2 the system shows the shortest (I), this is certainly a 
case of special interest. It becomes more curious when we notice that in one dimension there is not such an exception 
in the scaling analysis [10]. As shown in Fig. 3, at a = 3, (I') appears as a function of pLlnN, and the dynamic small 
world behavior is seen when this parameter exceeds 1. There is more discussion of these interesting points later. 

Figs.l, 2 and 3 can give us more information once we get an idea of L. For N large enough, we can use the following 
approximation, 



which gives 
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(12) 
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21n ^ _ln(W- 11 

(13) 

Some comments on this approximation: We have chosen the upper and lower limits for the integrals with some 
arbitrariness. For a ^ 3 and for N very large, it gives results good enough for the later discussion. Surely this 
approximation fails for a > 3, but in that region what is important is L stays finite as N goes to infinity. 
From Eqs. (9) to (13), we have for large enough N 

N 
ia=o-y, (14) 

L - - i(hW (15) 

1 — a 1 — 2 Q ~ 2 

L = 7 1 oa-i ^' (0 < « < 1. K a < 2) , (16) 

3 — a 1 — 2" 1 

L a = 2 -ln2^, (17) 

L^i V 3-«l_^./ 1 _2 a - 2 ) ^0, (2<a<3) (18) 

3 — a v y 



L Q=3 ^21nTV, (19) 

L -> finite, (a > 3) . (20) 

In the above equations, a = 0,2,3 come out as special points, a = 0, corresponds to the totally random network. 
Below a\ = D = 2, L^< a <2 is always proportional to TV. In Ref. [16], the authors proposed that for < a < a\ the 
system is in the random network phase, and a\ is the continuous phase transition point from the random network 
phase to the small-world phase. As a increases above cq, = D + 1 = 3, L is finite and independent of TV in the limit 
of TV — > oo. As a result, the system is virtually a regular network for a > cc c . 

Now we have the scaling relations shown in Figs.l, 2 and 3, and L as a function of TV given by Eqs. (14) to (20). 
In the following we shall discuss the system behavior with a starting from zero. 

(1) < a < a l c : As shown in Fig.l, when pL <C 1, (I') —> 0.5, indicating that the network is virtually regular. 
When pL increases beyond 1, the system begins to show a dynamic small world behavior, in the sense that 

<F)~(pLH~p-^V^, (21) 

or 

(j) ~ p-TJV 1 -T, (22) 

where 7 depends only on a. From the linear fit, at a — we have obtained that 7 « 0.309 and (I) ~ p-o.309^y0.69i^ 
Note that 7 is obtained from linear fit of limited data and cannot be exact, but our result, (3 = 1 — 7 « 0.691, is close 
to Kleinberg's result, /3 = 2/3. (We observe that in Ref. [11], (3 = 1/3.) As a increases, 7 increases, but (3 remains 
positive as a — > a \. With our present results, we are not able to give the full function of 7 (a), because near a = 2, 
(I) as a function of N severely deviates from a power law for relatively small N. 

(2) a = a l c = 2: At this point, when pL/\nN <C 1, (/') — ► 0.5, showing that it is a regular network. When 
pL/\nN increases above 1, exp I (I') x pL) turns into a linear function of pLj In TV, and it means 

(I 7 ) ~ In (pL/ In TV) /pL. (23) 

With L ~ ln2y-^ for TV — > 00, this gives 



or, in the limit of TV — > 00, 



— ■ In TV [In p + In (In 2)+ In TV -2 In (In 7V)1 , N 

(1) L^ !: ' v ^ 24 
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We have studied the case of < a < a\ , and the case of a > a\ will be studied below. We shall see that at a = a\ 
the smallest {I) is achieved. 

(3) a\ < a < o? c : Regularity dominates for pL <C 1. For pL > 1, dynamic small world effect arises. Similar to the 
case of < a < a\ , we have once again obtained 

Jl) ~ p-"f N 1 ^ , (26) 

and 7 tends to zero as a approaches a 2 . 

(4) a = c? c = 3: When pLlnN <C 1, (l 1 ) — ► 0.5, the network also shows dominating regularity. When pLlnN 
greatly exceeds 1, 

JFj-tpLlnNy" 1 . (27) 

Substitute L ~ 2 In TV into the above equation, then for large enough TV, we have 

"(O-p-TTV (In TVr 27 , (28) 

but the size of the networks used in the present study prevented us from getting a accurate estimate of 7. 

(5) a > a 2 c : Since in this case L stays finite when N — > 00, the system is believed to behave like a regular network. 
It is confirmed by the numerical calculation, which gives (I) nearly linear as TV. 



IV. SUMMARY AND DISCUSSION 

In this work, we investigate the navigation process on a variant of Watts-Strogatz (W-S) model embedded on a 
two-dimensional square lattice with periodic boundary condition. With probability p, each vertex sends out a long 
range link, and the probability that the other end of this link falls on a vertex at lattice distance r away decays as 
r~ a . Vertices on the network have knowledge of only their nearest neighbors. In a navigation process, messages are 
forwarded to a designated target, and the average actual path length (l) is obtained with varying a, p, and N. 

Our result is consistent with the existence of two phase transitions at or c = D = 2 (random network to small world 
network), and a 2 = D + 1 = 3 (small world network to regular network). For a < a 2 , and a ^ a\, it is found that 
(V) = (l)/N ~ f a (pL), where L is the average length of the additional long range links. This develops the scaling 
analysis in the works of Newman [15] and Zhu et al. [10]. Given pL > 1, dynamic small world effect is observed, and 
the behavior of f a at large enough pL gives (I) ~ p~T(") JV 1-7 ^'. When a = 0, 7 is close to 1/3, so (3 = 1 — 7 ~ 2/3, 
in agreement with Klcinbcrg's result of (/) oc N 2 / z for p = 1. As a increases, 7 increases (but stays below 1), and 
once a exceeds a*, 7 begins to decrease, and 7 approaches zero as a — > a 2 . At a = a\, this kind of scaling breaks 
down, and (V) can no longer be written as a function of pL. In this case we can still get (I) ~ In TV (hip + h\N) jp for 
N large enough. Note that only at a — a\, (I) grows as a polynomial of In A, and it is the closest point to the static 
small world effect [10]. At a — a 2 ., the scaling is (/') ~ / (pLlnN), and accordingly (I) ~ p~~<N (In TV) with large 
enough pLIn A. For a > a 2 , (I) is nearly linear with TV. 

It is reasonable that in social networks (like various other networks) the probability of connection falls as distance 
(in various senses, e.g., occupation) increases, and the apparently very small value of (/) in human society [1,2] suggests 
human society might have its exponent a being close to a\. 

The great success of the idea of small world has since motivated much effort in studying various dynamic processes 
based on the small world network model. The limited knowledge of the nodes of a network is an important limitation 
that has to be considered when studying navigation processes. Another interesting and important limitation is due 
to the fact that the links in a network are usually associated with "weights", as systematically studied in Refs. [17]. 
Further studies on link-weighted small-world model should help us gain insight in the navigation and other relevant 
phenomena in various artificial and natural networks, and help us design networks with higher efficiency. 
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Caption of figures 



Fig.l. (Color online) The reduced average actual path length (I') varies as pL for a = 0,1,2.5, where L is 
the average length of the additional long range links. The data collapse with each specific a contains curves with 
N = 200,400,600,800, 1000 respectively. On each curve with fixed N, p = 1.3 _i , where i = 0, 1,2, ...,47 (p has the 
same set of values on each curve in Fig. 2 and Fig. 3). The solid line y ~ x ~°- 309 is a guide to the eye. 

Fig. 2. (Color online) With given a = 2, plots show the relation between exp ((/')• pL J and pL at N = 

200,400,600,800,1000. 

Fig. 3. (Color online) The relationship between (I') and pL In N in log- log scale at a = 3 for N = 
200,400,600,800,1000. 
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